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Abstract 

We propose a modified Balitskii-Fadin-Kuraev-Lipatov equation from the 
viewpoint of the resummation technique, which satisfies the unitarity bound. 
The idea is to relax the strong rapidity ordering and to restrict phase space 
for real gluon emissions in the evaluation of the BFKL kernel. It is shown 
that the power-law rise of the gluon distribution function with the small 
Bjorken variable x turns into a logarithmic rise at x ^ 0. 
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It is known that the Bahtskii-Fadin-Kuraev-Lipatov (BFKL) equation 
[||] sums leading logarithms ln(l/x), x being the Bjorken variable, produced 
from the reggeon ladder diagrams with the rung gluons obeying the strong 
rapidity ordering. The gluon distribution function, governed by this equation, 
is found to increase at small x, which has been confirmed by the recent HERA 
data of the proton structure function F2 involved in deep inelastic scattering 
(DIS) 0. However, the increase is power-like, such that F2 and the DIS cross 
section cxtot rise as a power of x, i. e., as a power of s for x = Q"^ / s, Q being the 
momentum transfer and s the total energy. This behavior does not satisfy 
the Froissart bound atot < const, x In^ s, and violates unitarity. Hence, the 
BFKL equation can not be the final theory for small x physics. Though it has 
been expected that the inclusions of next-to- leading ln(l/x) and of higher- 
twist effects from the exchange of multiple ladders (pomerons) [|] may soften 
the BFKL rise, the attempts have not yet led to a concrete conclusion. In 
[^] a renormalization group (RG) improved double factorization formula was 
employed, and the obtained structure functions satisfy the unitarity bound. 
However, the logarithmic rise appears as the consequence of the RG, instead 
of BFKL, evolution. 

In this letter we shall propose a modified BFKL equation from the view- 
point of the resummation technique 0, which satisfies the unitarity bound. 
This technique was developed originally for the organization of double loga- 
rithms In^ Q, and applied recently to the all-order summations of various 
large logarithms contained in parton distribution functions. It has been 
shown that the known evolution equations, such as the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equation 0, which sums single logarithms 
InQ, the BFKL equation and the Ciafaloni-Catani-Fiorani-Marchesini 
(CCFM) equation |^, which embodies the above two equations, can be de- 
rived in a simple and unified way 

Following we have presented a straightforward application of the re- 
summation technique [0. Since the BFKL equation is independent of Q, its 
predictions are insensitive to the variation of Q. However, the HERA data 
for F2{x, Q"^) exhibit a stronger Q dependence 0: The rise at low x is slower 
for smaller Q. It is plausible that the experiments have not explored the 
region with x small enough to implement the BFKL evolution. To explain 
the data, the Q-dependent DGLAP equation is then combined in some way. 
For example, the CCFM equation was employed [Tl|]. Instead, we introduced 
the Q dependence into the BFKL equation by cutting off the longitudinal 
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component of the real gluon momentum at Q, when evaluating the BFKL 
kernel. The gluon distribution function derived from this Q-dependent BFKL 
equation leads to predictions of F2, which are well consistent with the data. 
Unfortunately, the predictions still rise as a power of x at x 0. Here 
we shall further modify the BFKL equation by relaxing the strong rapidity 
ordering in addition to restricting phase space of real gluon emissions. It is 
found that the resultant correction is negative, and the power-law rise of the 
gluon distribution function turns into a logarithmic rise at x ^ 0. 

First, we review the derivation of the BFKL equation using the resumma- 
tion technique @, [l^. The unintegrated gluon distribution function F{x, fcr), 
describing the probability of a parton carrying a longitudinal momentum 
fraction x and transverse momenta ky, is defined by 

F{x,kr) = - f ^f!^e-(^P"^-'^-^-) 

xWiP^a\F;i{y-,yT)F^^mP,^), (1) 

in the axial gauge n ■ A = 0, n being a gauge vector with 7^ 0. The ket 
\p, a) denotes the incoming proton with the light-like momentum p'^ = p'^5^^ 
and spin a. An average over color is understood. F^ is the field tensor. 

For a fixed parton momentum fc"*", F varies with p"*" implicitly through 
X = k^/p'^, and we have p'^dF/dp'^ = —xdF/dx. That is, the derivative of 
F with respect to x is related to the derivative with respect to p^. Because 
of the scale invariance in n of the gluon propagator —iN^'^{l)/P with 

iV^-=^/^- + ^2^^^ 2 

n ■ I [n- ly 

the logarithmic corrections to F must depend on the ratio {k ■ ny jr? cx 
(p ■ n^jr?^ if assuming n = (n"*", ra", 0). Hence, there exists a chain rule 
relating p^d/dp^ to d/dn, 

^ dp~^ V ■ n ^ drip ' 

V(3 = (5/3+ being a vector along p. The operator d/drip applies to a gluon 
propagator, giving 

J^N"""' = ^—irN^"' + F'N'''^) . (4) 

drip n ■ I 
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The loop momentum l'^ (l'^') contracts with a vertex in F, which is then 
replaced by a special vertex -U/j = n^vp/iv -nn-l). This special vertex can be 
read off the combination of Eqs. (0) and (^). 

The contraction of l'^ leads to the Ward identity, 



-im^<{k) -iN'^^'ik + l) 



F {k + iy 



(5) 



^fiux being the triple-gluon vertex. Summing all the diagrams with different 
differentiated gluons, those embedding the special vertices cancel by pairs, 
leaving the one in which the special vertex moves to the outer end of the 
parton line We then obtain the derivative, 

- x-^F{x, kr) = 2F{x, kr) , (6) 

described by Fig. 1(a), where the new function F contains one special vertex 
represented by a square. The coefficient 2 comes from the equality of the new 
functions with the special vertex on either of the two parton lines. Equation 
(ID is an exact consequence of the Ward identity without approximation [^]. 
An approximation will be introduced, when F is related to F by factorizing 
out the subdiagram containing the special vertex, such that Eq. (|]) reduces 
to a differential equation of F. 

It is known that factorization holds only in leading regions. The leading 
regions of the loop momentum / flowing through the special vertex are soft 
and hard, since the factor 1/n ■ I with 7^ in Vj3 suppresses collinear 
enhancements 0. For soft and hard /, the subdiagram is factorized according 
to Figs. 1(b) and 1(c), respectively. A straightforward evaluation shows 
that the contribution from the first diagram of Fig. 1(c) is less important 
in the considered region with small k~^. Therefore, we drop Fig. 1(c), and 
concentrate on Fig. 1(b). The color factor is extracted from the relation 
fabcfbdc = —NcSad, whcre the indices a,b, . . . have been indicated in Fig. 1(b), 
and Nr = 3 is the number of colors. The factorization formula is written as 



soft 



X 



2m5{f)F{x + |kT + ItI) + ^(x, kr) 



(7) 
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where iNc comes from the product of the overall coefficient —i in Eq. (j^) and 
the color factor —N^ extracted above. The triple-gluon vertex for vanishing 
/ is given by 

r^i^A = ~gf,uxpx - guxxp^, + 2gx^,xp^ . (8) 

The denominator —2xp ■ I is the consequence of the eikonal approximation 
for the gluon propagator, {xp — l)'^^ —2xp ■ I. The first term in the brackets 
corresponds to the real gluon emission, where F{x + l~^/p'^, IkT + lrl) implies 
that the parton coming out of the proton carries the momentum components 
xp'^ + and k-r + Ir in order to radiate a real gluon of momentum /. The 
second term corresponds to the virtual gluon emission, where the 6 function 
sets the upper bound of It to to ensure a soft momentum fiow. 

It can be shown that the contraction of p with a vertex in the quark 
box diagram the partons attach, or with a vertex in the gluon distribution 
function, leads to a contribution down by a power 1/s, s = {p+qY, compared 
to the contribution from the contraction with vp. Following this observation, 
Eq. (|^) is reexpressed as 



soft 



2 72 \ 



+^-^r^Fix,kT) 



(9) 



The eikonal vertex comes from the last term xpi, (divided by xp'^) in 
Eq. (j^), and the eikonal propagator 1/v-l from l/{xp-l), which is represented 
by a double line in Fig. 1(b). The remaining metric tensor g^^ has been 
absorbed into F. 

To derive the conventional BFKL equation, we simply assume the strong 
rapidity ordering, x+l^ /p'^ ^ x, for the real gluon emission, namely, approx- 
imate F{x+l'^ Ip'^ , \kT+'^T\) by its dominant value F{x, \k.T+'^T\)- Performing 
the integrations over l~ and Z+ to infinity, and substituting F ^ Fsoit into 
Eq. (P), we arrive at 



dF{x, kr) _ f d% 

a. 



d\n{l/x) "J T^ll 



F{x,\\iT + M)-0{kl-ll)F{x,kT) , (10) 



with as = Ncas/T^i which is the BFKL equation. To obtain the Q-dependent 
BFKL equation, we truncate Z"*" at Q/\/2 for the real gluon emission in 
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Eq. (I), giving 



soft 



as_ 

y 





/ 


I TT 


Jo 



dl 



21+ 



n 



xFix + t/p+, |kT + ItI) - ^-^^^^Fix, kT) 



.(11) 



The motivation is that the vanishing of F{x + 1+ /p+) at large momentum 
fraction constrains 1+ to go to infinity. Applying the strong rapidity ordering, 
Eq. (0) leads to 



dF{x, kx) 
ciln(l/a;) 



d^T 



F{x, \kT + ItI) - ^(4 - ^T)^(a;, kr] 



-a. 



dHTF{x,\kT + lT\) 



TT 



(12) 



for the choice n = (1, 1, 0). The last term, corresponding to the upper bound 
of l+, can be regarded as bringing in some higher-power contributions, which 
moderate the BFKL rise at low Q. 

As stated before, the Q-dependent BFKL equation (|12|), though phe- 
nomenologically successful, predicts a power-law rise for the gluon distri- 
bution function at x — >■ 0, which violates unitarity. We point out that the 
assumption of the strong rapidity ordering may be the cause for the unitarity 
violation. For most values of 1+ , F{x + l+/p+, Ikr + lrl) is much smaller than 
F{x, \]<.t + It\), and replacing the former by the latter in the whole integration 
range of 1+ overestimates the contribution from real gluon emissions. Hence, 
we shall not adopt the assumption and employ Eq. 
variable change = yy/xp+, we obtain 



11) directly. Using the 



dF{x, kj) 
d\n.(l/x) 



rdHr 


IL 


1 TT 





9{kl- f 



^'F{x,kT) , (13) 



where the upper bound of y is determined by the kinematic relation Q ^ 
\/2xp+. Equation (0) is the modified BFKL equation we shall investigate 
in more details below. 
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An initial condition F{xo, kx) = F^'^^xq, kx) must be assumed when solv- 
ing Eq. ([T3|), Xo being the initial momentum fraction below which F evolves 
according to the BFKL equation. For instance, a "flat" gluon distribution 
function ITTII 



F^^\x,kj 



1 GeV' 



x] 



'exp[-4/(l GeV^)] 



(14) 



with Ng a normalization constant, has been proposed. Hence, the initial 
function F^^\x + yy/x, Ikr + lrl) should be substituted for F{x + yy/x, \k.T + 
\t\) in Eq. (p!3D as x + yy/x > xq- 

In order to simplify the analysis, the 6 function for the virtual gluon 
emission is replaced by 9{Ql — Z^) [|10|, where the parameter Qq can be 



determined from data fitting. This modification is acceptable, because the 
virtual gluon contribution only plays the role of a soft regulator for the real 
gluon emission, and setting the cutoff of It to Qq serves the same purpose. 
Reexpressing the integrand F{x + yy/x, |kr + 1^1) as 



F{x, \kT + 1t|) + \f{x + yV^, \kT + ItI) - F{x, |kr + It|) 



(15) 



and Fourier transforming Eq. 
dF{x,b) 



into the b space conjugate to kx, we derive 



dln{l/x) 



S{b,Q)Fix,b) 



+2as{l/b)Qb / dyKi{yQb)[F{x + y^,b)-F{x,b)] 



(16) 



where 



Sib, Q) = -2a,(l/6) [ln(Qo&) + 7 - In 2 + K^iQb)] (17) 

comes from the combination of the first term in Eq. (|15|) and the virtual 
gluon emission term. Kq and Ki are the Bessel functions, and 7 the Euler 
constant. The argument of ag has been set to the natural scale 1/6. The 



second row of Eq. (IC) is identified as the negative correction from relaxing 
the strong rapidity ordering. 

Before solving Eq. (|1^), we extract the behavior of F analytically. Sub- 
stituting a guess F oc x^^ into Eq. (0), A being a parameter, we obtain 

1 r / , ^\ -A 

X = S + 2a,Qb [ dyK^iyQb) 
Jo 



' X + y^/x^ 



X 



(18) 
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It can be numerically verified that a solution of A, < A < exists for 
X < xq. That is, F increases as a power of x, consistent with the results 
from the conventional and Q-dependent BFKL equations. The correction 
term diverer 0, and no solution of A is allowed, implying that F can 

not maintain the power-law rise at extremely small x. We then try another 
guess F\x^Q oc ln(l/a;) with a milder rise. It is found that the correction 
term, increasing as ln^'^(l/a;), only slightly dominates over the first term 
SF oc ln(l/x) (the derivative term —xdF/dx oc 1 is negligible), and Eq. (|1^) 
holds approximately. At last, we assume -F|x-*o oc const, as a test. In 
this trial the first term becomes dominant, and the correction term and the 
derivative term vanish, i.e., no const. 7^ exists. These observations indicate 
that F should rise as ln(l/a:) at most when x approaches zero. We conclude 
that the modified BFKL equation predicts a rapid power-like rise of F for 
X < Xq and a milder logarithmic rise at x ^ 0. 

The behaviors of other quantities can be deduced. The gluon density, 
defined by 

/•Q fl^hrp roo _ 

xg{x,Q^)= -F{x,kT)=2Q db.JiiQb)F{x,b) , (19) 

Jo 71 Jo 

possesses a similar dependence on x, where Ji is the Bessel function. The 
structure function F2 is written, in terms of the /cT-factorization theorem, as 

F2{x, Q') = f'^ f ^H{x/^, kr, Q)F(e, kr) , (20) 

Jx t, J TV 

where the hard scattering subamplitude H denotes the contribution from 
quark box diagrams. Because of F{^,kT) < ln(l/^) at ^ — > 0, F2 rises as 
ln^(l/x) at most as indicated by the ^ integration, and thus satisfies the 
unitarity bound. 

We present the explicit numerical results of the gluon density xg in Fig. 2 
with Ng ^ 3.65 and Qo ~ 0.3 |jlO[, which, obviously, behave in the way 
stated above. The deviation of the curve from a straight line at x — > is 
the consequence of relaxing the strong rapidity ordering. The slower rise 
of the curve for a smaller Q is due to the truncation of /"^ at Q/\/2. The 
comparision of our predictions for F2 with the experimental data and the 
connection of the modified BFKL equation ([T3| ) to the other approaches in 
the literature will be published elsewhere. 

This work was supported by the National Science Council of R.O.C. under 
the Grant No. NSC-87-2112-M-006-018. 
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Figure Caption 



FIG. 1. (a) The derivative —xdF/dx in the axial gauge, (b) The soft struc- 
ture and (c) the ultraviolet structure of the 0{as) subdiagram containing the 
special vertex. 

FIG. 2. The dependence of xg on x for Q"^ = 8.5 GeV^ (dotted hne), 15 
GeV^ (dashed line), and 20 GeV^ (sohd line). 
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